Trapping and transmission of matter-wave solitons in a collisionally inhomogeneous 

environment 
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We investigate bright matter-wave solitons in the presence of a spatially varying scattering length. 
It is demonstrated that, even in the absence of any external trapping potential, a soliton can be 
confined due to the inhomogeneous collisional interactions. Moreover we observe the enhanced 
transmission of matter-wave solitons through potential barriers for suitably chosen spatial variations 
of the scattering length. The results indicate that the manipulation of atomic interactions can 
become a versatile tool to control the dynamics of matter waves. 
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Introduction. The recent developments in the field 
of atomic Bose-Einstein condensates (BECs) |l| has in- 
spired, among others, many studies on their nonlinear 
excitations. Especially, as far as matter-wave solitons are 
concerned, dark |2|, bright and gap |f| solitons have 
been observed experimentally and studied theoretically. 
Atom optical devices such as the atom chip [fj offer the 
possibility to control and manipulate matter-wave soli- 
tons. Their formal similarities with optical solitons indi- 
cate that they may be used in applications similarly to 
their optical counterparts Q. 

Typically bright (dark) matter-wave solitons are 
formed in BECs with attractive (repulsive) interatomic 
interactions, i.e., for atomic species with negative (pos- 
itive) scattering length a. Employing magnetically- 
induced Feshbach resonances both the magnitude and 
sign of the scattering length can be changed by tuning 
the external magnetic field (see e.g. || and also 0, 0] 
where the Feshbach resonance in Li condensates was 
used for the formation of bright matter- wave solitons). 
These studies paved the way for important experimen- 
tal discoveries, such as the formation of molecular BECs 
H and the revelation of the BEC-BCS crossover [Icj . 
From the theoretical viewpoint, it was predicted that a 
time-dependent modulation of the scattering length can 
be used to prevent collapse in higher-dimensional attrac- 
tive BECs pj| , or to create robust matter- wave breathers 
[lij l . Adding to a constant bias magnetic field a gradient 
in the vicinity of a Feshbach resonance allows for a spa- 
tial variation of the scattering length over the ensemble 
of cold atoms thereby yielding a collisionally inhomoge- 
neous condensate. Due to the availability of magnetic 
and optical (laser-) fields the external trapping poten- 
tial and the spatial variation of the scattering length can 
be adjusted independently. Moreover hyperfine species 
with the magnetic quantum number Mjr = do not 
feel a potential due to the magnetic field but experience 
magnetically-induced Feshbach resonances (see e.g. ) ■ 
In this case the external potential is formed exclusively 
by an optical dipole potential and the magnetic field con- 



figuration is responsible for the spatially dependent scat- 
tering length. Recently this has been exploited to study 
the properties of cold atomic gases in a collisionally in- 
homogeneous environment (CIE) [13, H^ . 

Here we investigate the dynamics of bright matter- 
wave solitons of a quasi one-dimensional (ID) 7 Li BEC 
[jj in a CIE demonstrating the appearance of unex- 
pected phenomena which make the spatial manipulation 
of the scattering length a versatile tool. Firstly we show 
that the inhomogeneity of a induces an effective confin- 
ing potential felt by a bright matter-wave soliton, even 
in the absence of any external trapping potential. The 
proposed scheme leads to a collision-induced breathing 
soliton, which oscillates due to the effective confinement 
and is periodically compressed when passing through the 
region of a large scattering length. As a second prototyp- 
ical situation for a CIE, we consider the transmission of 
a soliton through a potential barrier underneath which a 
suitably chosen spatially varying scattering length a(x) 
is present. This setup allows to enhance the transmission 
of the soliton, i.e., the barrier becomes more transparent 
compared to the case of a spatially independent atom- 
atom interaction. More specifically, we use the partic- 
ular form a(B) for a 7 Li condensate near a correspond- 
ing Feshbach resonance (see below). This species has al- 
ready been used in order to prepare bright matter-wave 
solitons. We believe that the collision-induced trapping 
and transmission of matter-wave solitons reported here, 
are two generic phenomena illustrating that collisionally 
inhomogeneous matter-waves exhibit a number of inter- 
esting and fundamentally new features that could also be 
relevant for future applications. 

Collision-induced trapping. To specify our setup we 
choose the magnetic field dependence of the scattering 
length a of a 7 Li condensate provided in Ref. Q (see also 
Ref. Q). However, we emphasize that this case serves 
only as a typical example of such a curve, used for con- 
creteness. Let us focus on the regime < B < 590 G, 
which is far left from the Feshbach resonance at 720 G, 
where the scattering length is small and the inelastic col- 
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FIG. 1: The spatial variation of the scattering length for B = 
450 + ex (G) and magnetic field gradients e = (dotted line) , 
e = 5 G/fj,m (solid line) and e — 25 G/fj,m (dashed line). 



lisional loss of atoms is practically negligible. We have 
a(B) < for 150 G < B < 520 G and a(B) > else- 
where. At B w 352 G, the scattering length reaches a 
minimum with a m —0.23 nm. The experimentally ob- 
served quasi- ID bright matter- wave solitons 0, 0] have 
been created in the above-given regime of negative scat- 
tering length. 

Let us assume a magnetic field configuration B = 
Bq + ex (G), where e is the field gradient and Bq = 450 G, 
being far from the position of the resonance. The gradi- 
ent chosen is of the order of a few tens of G//im which can 
be experimentally realized in microscopic matter-wave 
devices such as the atom-chip 6J. Other resonances and 
species might require much smaller values of the gradient 
to implement a significant change of the scattering length 
on the scale of the condensate. To extract the spatial de- 
pendence of the scattering length we apply a fifth-order 
polynomial fit to the data for a{B) given in Ref. Q. 
The result is shown in Fig. 1 for two nonzero values of 
e. The function a(x) possesses a single minimum and the 
dependence on x is, as expected, much more dramatic for 
a larger value of e. For e = we have a — —0.182 nm. 

The evolution of an untrapped, quasi-lD bright matter- 
wave soliton in a CIE is described by the following nor- 
malized Gross-Pitaevskii (GP) equation: 
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d 2 x ^ - g(xM 2 ip. 



(1) 



Here, tp is the mean field wavefunction (with the den- 
sity | ip\ 2 measured in units of the peak density no), x 
is given in units of the healing length £ = h/^/nogom 
(where go — 2hioj_a(Bo) and luj_ is the confining fre- 
quency in the transverse direction), and the time unit 
is £/c (where c = y/nogo/m is the Bogoliubov speed of 
sound). Finally, the spatially dependent nonlinearity is 
given by g(x) — a(x)/a(Bo), with a{x) = a(Bo + e£x) 
(note g{x = 0) = g(e = 0) = 1). Typically, for a quasi- 
lD 7 Li condensate with lj± — 2ir x 1000Hz and no = 10 9 
m _1 , the healing length and speed of sound amount to 




FIG. 2: (Color online) Spatio-temporal contour plot of the 
density of a bright matter wave soliton for Bo = 450G 
and magnetic field gradients e = 5G/fim (left panel) and 
e = 25G//xm (right panel). The dashed line in the left panel 
corresponds to the analytical prediction according to Eq. |QJ|. 



£ = 2«m and c = 4.6mm/s, respectively. 

Introducing the transformation ip — u/yfg in the re- 
gion g(x) > we reduce Eq. JIJ to the following per- 
turbed nonlinear Schrodinger (NLS) equation, 
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\d 2 x i 
2 x 



\u\ 2 u 



R(u), 



(2) 



with a perturbation R(u) = -4- \n(y/g)d x u + O (L| / L^) . 
The first term of R is of the order O (Ls/Lb), where 
Ls and Lb are the characteristic spatial scales of the 
soliton and of the inhomogeneity due to the magnetic 
field gradient, respectively. In the case R = (<J4> e = 0), 
Eq. (J2J has a bright soliton solution of the form , 



u(x, t) — ?7sech[?7(a; — xq)] exp[i(fciE — <f>(t)], 



(3) 



where r\ is the amplitude and inverse width of the soliton, 
and xq is the soliton center. The parameter k — dxo/dt 
defines both the soliton velocity and wavenumber, and 
4>(t) = (l/2)(k 2 — rj 2 )t is the soliton phase. For the above 
mentioned typical values of the parameters, and for n = 1 
(a soliton with 10 3 atoms and width 2£ = 4/im), it is clear 
that Ls/Lb = 2£e/£?o. This means that for sufficiently 
small values of the magnetic field gradient, e.g., for e = 
5 G/um, the perturbation R is of order O(10 -2 ). In 
such a case, we may employ the adiabatic perturbation 
theory for solitons [ljj, to obtain the following equation 
of motion for the soliton center, 

d 2 x 



dt 2 



dVeB 

dx ' 



V e s(x ) 



(4) 



With g{x) being proportional to a{x), it is obvious 
(see Fig. that the soliton "feels" a collision-induced 
effective confinement potential V e s_, although there is no 
external trapping potential [see Eq. IjTjl]. It is then natu- 
ral to expect that the collision-induced confinement leads 
to oscillations of the soliton if it is displaced from the 
minimum of the confinement potential. This has been 
verified by direct numerical integration of the GP Eq. 

for e = 5 G//im and initial condition ip = sech(a;) 
(i.e., the soliton is initially at xq(0) = where B = Bo). 
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The result is shown in the left panel of Fig. 2, where 
the spatio-temporal contour plot of the soliton density 
is directly compared to the analytical prediction of Eq. 
Q (dashed line); the agreement between the two is ex- 
cellent. We observe that the matter- wave is periodically 
compressed whenever it reaches the region of large scat- 
tering lengths (x ss — 18/im with a w — 0.23nm), thus ex- 
hibiting a pronounced breathing behavior. This results 
in a spontaneous and robust breathing behavior of the 
matter-wave soliton in the collisionally inhomogeneous 
environment. 

For significantly larger field gradients, e.g., e = 25 
G//im, the soliton width becomes comparable to the 
magnetic length scale L b ■ This means that the perturba- 
tion R in Eq. J5J is now of the order L S /L B = O(10 _1 ) 
and nonadiabatic effects are expected to be significant. 
The numerical integration of the GPE confirms this ex- 
pectation and reveals that although the soliton is still 
confined and performs corresponding oscillations, its evo- 
lution is nonadiabatic, i.e., emission of small amplitude 
wave radiation is observed. The corresponding results are 
shown in the right panel of Fig. 2: Larger field gradients 
lead to oscillations with a higher frequency and smaller 
amplitude. 

Collision-induced transmission. In our second setup 
we consider the scattering of a bright matter- wave soliton 
(see Eq. (jSJ for t = and rj = 1) of an external potential 
barrier thereby comparing the results for homogeneous 
and inhomogeneous atomic interactions. An important 
quantity in this context is the transmission coefficient T. 

In order to compare the transmission in the above 
mentioned cases, the incoming and outgoing scattering 
environments should be identical. This means that the 
scattering length should asymptotically, i.e. outside the 
range of the barrier, take on the same values. We guar- 
antee this by employing a localized inhomogeneous mag- 
netic field of the form, 

B(x) = - [(Bi + B 2 ) + (B 1 - B 2 ) t&nh(w(x - x B ))} , (5) 

where the parameters w and xb characterize the inverse 
width and location of the region of inhomogeneity, re- 
spectively, while the field values B\ and B 2 are chosen 
to obtain equal scattering lengths a(B 1 ) = a(B 2 ) suffi- 
ciently far from the barrier. This value will also be used 
for the homogeneous case. In the case of 7 Li it is straight- 
forward to find the values B\ and B 2 , due to the conve- 
nient form of the function a{B) in the considered range 
of field strengths (see also Fig.l). Here we use B\ = 450 
G and B 2 = 265 G for which a = -0.182 nm (other 
choices are, of course, equally possible and lead to sim- 
ilar results). We finally note that the above mentioned 
field configuration can be realized by a multi-wire setup 
or a current density flowing in a half plane augmented by 
a homogeneous bias field. 

Our potential barrier is assumed to be of the form 




FIG. 3: (Color online) The scattering of a soliton initially 
placed at x = off a barrier (shaded area) located at xb = 
— 20/^m is shown; the transmitted and reflected parts of the 
soliton (accordingly labeled) are shown for t — 6ms for the 
inhomogeneous case. The spatial dependence of the scattering 
length is also shown. Note that for the homogeneous case 
a = -0.182nm. 

Vb[x) — Vosech 2 (a(x — xb)), where Vo, a -1 and xb are 
the barrier's amplitude, width and location, respectively 
(note that the inhomogeneity is centered at the same po- 
sition where the barrier is located, i.e., at x = Xb). In 
the following we assume Vo = 1 and a' 1 — 1/2, i.e., the 
width of the barrier is half the soliton width, so as to 
avoid the classical Ehrenfest regime. 

The setup is illustrated in Fig. 3: The initial (t = 0) 
form of the soliton, as well as its transmitted and re- 
flected parts (for the inhomogeneous case at t = 6ms) 
are depicted and labeled respectively. At the location 
of the barrier (shaded region), and for the collisionally 
inhomogeneous case, there exists a local spatial change 
of the scattering length, which is obtained by inserting 
the magnetic field in Eq. © into the function a(B) (uti- 
lizing the fifth-order polynomial fit used in the previous 
section). The spatial dependence of the scattering length 
is also shown in Fig. 3; as discussed above, the scattering 
length takes on equal values far from the barrier for both 
the homogeneous and inhomogeneous case. 

We have numerically integrated the GPE incorporat- 
ing the potential barrier term Vb(x)ip(x,t) to determine 
the transmission coefficient T. The results presenting 
T as a function of the width of the inhomogeneity (left 
panel) and the soliton's incident velocity (right panel) 
are shown in Fig. Generally, for a fixed soliton ve- 
locity, or width of the inhomogeneity, the transmission 
T in the inhomogeneous case is always larger than the 
one in the homogeneous case. Particularly, for k = 1.1c 
and wt; = 1.3, the relative difference of the transmis- 
sion T for the two cases becomes maximal, being ss 15%. 
This result clearly demonstrates that the transmission of 
a matter-wave soliton may be enhanced in the presence 
of a spatially dependent collisional interaction. 

Conclusions. We have explored the dynamics of bright 
matter-wave solitons subject to a spatially varying non- 
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FIG. 4: The transmission T as a function of the inverse width 
w of the inhomogeneity, in units of £~ (left panel), or the 
soliton's incident velocity k, in units of c (right panel) . Dashed 
and solid lines correspond to the collisionally homogeneous 
and inhomogeneous cases, respectively. 



linearity, which can be realized by means of an external 
inhomogeneous magnetic field on top of a bias field in or- 
der to be close to a region of strong changes of the scatter- 
ing length, such as a Feshbach resonance. It was demon- 
strated that a confinement (trapping) of the matter- wave 
can be achieved solely on basis of the spatial change of 
the collisional interaction, i.e., without the presence of 
an external trapping potential thereby creating a breath- 
ing matter-wave state in the collisionally inhomogeneous 
environment. In an adiabatic regime, such a state could 
be well described within the realm of soliton perturba- 
tion theory, while for abrupt variations of the scattering 
length, radiative emissions are non-trivial resulting in en- 
ergy losses and hence shorter period oscillations. Using 
a localized spatial variation of the scattering length, we 
have shown that the transmission of matter-wave solitons 
through a barrier can be enhanced by suitably manipu- 
lating the collisional properties of the condensate in the 
vicinity of the potential barrier. Collisionally inhomoge- 
neous environments therefore hold considerable promise 
in the effort to control and manipulate matter-waves in 
experimental applications. Interesting future directions 
may involve combining this type of effective potential 
with optical potentials in order to examine transmission 
properties and symmetry breaking or nonlinear trapping 
phenomena similarly to the recent works of [lij . An in- 
teresting variation on these themes that the present set- 
ting offers is, among others, that the effective potential 
is asymmetric and can hence lead to a modified dynami- 
cal picture [19| in comparison to the symmetric potential 
case [l^ . Such studies are currently in progress and will 
be reported in future publications. 
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